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ABSTO^T 

The  science  of  navigation  has  been  augmented  by  the  use  of  radio  aids,  of  which 
a  recently  considered  one  has  been  that  of  utilizing  radio -frequency  energy  reflected 
from  the  moon.  Two  characteristics  of  such  a  signal  have  been  examined  to  deter¬ 
mine  their  applicability  in  a  terrestrial  navigation  system;  first,  the  Doppler  fre¬ 
quency  shift  undergone  by  the  signal,  and  second,  the  time  rate  of  change  of 
Doppler  shift. 

The  ratio  of  Doppler  frequency  shift  to  the  time  rate  of  change  of  Doppler  shift 
is  a  function  of  the  longitude  of  the  observer,  and  the  ratio  of  Doppler  shift  to  the 
sine  of  the  local  hour  angle  of  the  moon  is  proportional  to  the  latitude  of  the 
observer. 

The  probable  error  in  position  on  the  surface  of  the  earth,  resulting  from  the 
probable  errors  in  the  parameters  of  an  idealized  system,  range  from  a  few  meters 
at  latitudes  near  90  degrees  to  approximately  100  kilometers  at  equatorial  latitudes; 
however,  nagivation  to  within  line  of  sight  (10  km)  of  a  given  position  maybe 
accomplished  over  70  percent  of  the  hemisphere  illuminated  by  the  moon  at  a  given 
time.  The  position  error  is  a  function  of  position  as  well  as  time  and  is  relatively 
independent  of  the  transmitter  frequency  employed.  The  ultimate  limitations  in 
position  error  are  those  due  to  the  fundamental  uncertainty  in  relative  velocity 
between  earth  site  and  moon  center  and  to  the  probable  error  in  transmitter  fre¬ 
quency.  The  ratio  of  probable  error  in  transmitter  frequency  to  the  transmitter 
frequency  should  be  of  the  order  of  1  part  in  10^‘  if  the  full  capabilities  of  the  sys¬ 
tem  are  to  be  realized. 

Two  experiments  were  performed  at  field  sites  of  the  Laboratory  as  part  of 
this  investigation  and  the  results  indicate  that  the  position  errors  described  here 
are  of  the  right  order  of  magnitude. 


PROBLEM  STATUS 

This  report  completes  one  phase  of  the  problem;  work  on  other  phases 
continues. 


AUTHORIZATION 

NRL  Problem  R06-13 
Projects  NR  417-000  and  NR  417-001 


Manuscript  submitted  August  7,  1957 


MOON  DOPPLER  NAVIGATION 
[Secret  Title] 


INTRODUCTION 


The  first  successful  radar  contact  with  the  moon  was  reported  in  1946  by  the  United 
States  Army  (1).  The  years  following  this  event  have  been  marked  by  increasing  interest 
in  the  reflection  of  radio  waves  from  the  moon.  It  has  been  proposed  that  techniques 
developed  in  the  course  of  experiments  conducted  in  this  branch  of  radio  astronomy  be 
employed  to  achieve  the  following  ends: 

1.  determine  the  earth -moon  distance, 

2.  measure  the  relative  velocity  of  earth  and  moon, 

3.  obtain  the  velocity  of  propagation  of  radio  waves  in  the  earth-moon 
space,  and 

4.  investigate  the  possible  applications  of  data  obtained  during  experiments 
of  this  nature. 


This  report  had  its  origin  in  4.  Its  purposes  are  to: 

1.  study  the  Doppler  shift  in  the  received  frequency  of  a  radio-frequency 
transmission  from  earth  to  moon  to  earth  as  a  function  of  earth  rotation 
and  orbital  motion  of  the  moon,  and 

2.  investigate  the  possibility  of  establishing  a  terrestrial  navigation  system 
based  on  contours  of  constant  Doppler  frequency. 

The  navigation  system  to  be  analytically  studied  is  one  where  a  single  site  is  used 
for  transmission  and  reception.  This  arrangement  is  termed  “monostatic.  ”  Following 
the  analytic  treatment,  limited  experimental  results  are  given  for  a  bistatic  arrangement 
where  the  transmitting  and  receiving  sites  are  separated  by  25  miles. 

The  term  “navigation”  implies  the  calculation  of  position  and  direction  with  minimum 
error.  If  there  is  a  given  error  in  the  location  of  a  fixed  site  on  the  surface  of  the  earth, 
and  the  Doppler  measuring  system  is  free  of  error,  the  uncertainty  in  the  location  of  the 
position  of  the  Doppler  measuring  site  relative  to  the  position  of  the  fixed  site  is  equal  to 
the  given  error.  If,  however,  the  Doppler  measurement  introduces  error,  the  uncertainty 
in  location  is  amplified.  The  errors  introduced  by  the  Doppler  measuring  technique,  and 
the  magnitude  of  the  resulting  amplification  of  position  uncertainty  on  the  surface  of  the 
earth,  will  be  examined. 
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THEORY 

When  a  radio -frequency  signal  is  transmitted  from  the  earth,  reflected  by  the  moon  and 
received  on  the  earth,  it  undergoes  a  double  Doppler  frequency  shift,  one  on  going  and 
another  on  the  return  journey.  The  total  shift  (Appendix  B)  is  given  by  the  expression 


Af  =  f,  - 

(1) 

"  * (¥) 

(2) 

■where 


f  is  the  frequency  received  on  the  earth 
is  the  frequency  transmitted  from  the  earth 
V  is  the  relative  velocity  of  the  site  on  the  earth  and  the  moon 
c  is  the  velocity  of  propagation  in  the  earth -moon  space 
Af  is  the  Doppler  frequency  shift. 

The  Doppler  frequency  shift  is  a  consequence  of  the  earth’s  rotation  and  the  moon  s  orbital 
motion.  It  is,  as  shown  in  Eq.  (2),  proportional  to  the  relative  velocity  between  the  moon 
and  the  receiving  site  on  the  earth. 

A  view  of  the  earth-moon  system,  from  a  point  above  the  earth’s  north  geographic 
pole,  is  shown  in  Fig.  1, 

where 


H  is  the  earth’s  north  geographic  pole 
A  is  a  receiving  site  located  on  the  earth’s  equator 
L  is  the  moon 
R  is  the  earth’s  radius 

0}  is  the  earth’s  angular  velocity  of  rotation 

LHA  is  the  angle  between  HA  and  the  line  connecting  earth  and  moon 
centers  (this  is  called  the  local  hour  angle  of  the  moon). 


If  the  moon  lies  in  the  earth’s  equatorial  plane  (declination  0  degrees),  the  relative 
velocity  of  the  site  A  with  respect  to  the  moon  (due  to  the  rotation  of  the  earth)  is 
-O)  R  sin  LHA  (a  positive  velocity  exists  when  the  distance  between  earth  site  and  moon 
center  is  diminishing).  Let  a  plane  be  constructed  so  that  it  includes  the  earth  s  polar 
axis  and  the  line  connecting  earth  and  moon  centers.  All  points  on  the  earth  lying  in  a 
plane  drawn  through  site  A,  parallel  to  the  plane  just  described,  will  have  equal  relative 
velocities  with  respect  to  the  moon  (since  the  earth  is  a  rigid  body).  However,  this  is 
strictly  true  only  if  the  moon  is  stationary,  relative  to  the  earth’s  center  and  the  fixed 
stars  and  at  an  infinite  distance.  This  condition  is  assumed  to  avoid  introducing  complex 
factors  at  this  point.  Figure  2  is  a  view  of  the  earth  from  a  direction  mutually  perpen¬ 
dicular  to  the  earth’s  polar  axis  and  the  line  connecting  earth  and  moon  centers,  and 
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\  TO  MOON 

-cjR  sin  LHA 


5.  1  -  Relative  velocity  of  site  A 
with  respect  to  moon 


NORTH  GEOGRAPHIC  POLE 


EQUATOR 


SOUTH  GEOGRAPHIC  POLE 


Fig.  2  -  Semicircles  of  equcil  velocity  on  the  earth’s 
surface  with  respect  to  moon,  viewed  from  zenith  of 
geographical  point  having  latitude  0  degree  and  local 
hour  angle  270  degrees 
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illustrates  the  intersection  of  planes  of  equal  velocities  and  the  earth’s  surface.  In  this 
view,  where  the  meaningful  curves  of  equal  velocity  are  coaxial  semicircles,  the  point  of 
maximum  velocity  is  the  point  of  intersection  of  the  earth’s  surface  and  the  line  through 
the  centers  of  the  semicircles,  that  is,  the  point  on  the  earth’s  surface  having  latitude 
0  degree  and  local  hour  angle  270  degrees.  Proceeding  around  the  earth  from  this  point 
one  encounters  semicircles  of  diminishing  velocity  until  the  meridian,  whose  plane  includes 
the  moon,  is  reached  where  the  relative  velocity  with  respect  to  the  moon  is  zero.  Since 
the  Doppler  frequency  shift  is  proportional  to  the  relative  velocity, these  curves  are  con¬ 
tours  of  constant  Doppler  shift.  They  may  be  thought  of  as  a  fixed  grid  or  cage  of  constant 
Doppler  frequency  lines  in  which  the  earth  rotates. 

To  introduce  the  concept  of  moon  Doppler  navigation.  Fig.  3  is  presented.  This  is  a 
view  of  the  earth -moon  system  from  a  point  above  the  earth’s  north  geographic  pole.  The 
parallel  lines  represent  the  intersection  of  equal  velocity  planes  with  the  surface  of  the 
earth.  The  relative  velocity,  with  respect  to  the  moon,  of  an  observer  moving  from  point  A 
to  B  to  C  will  vary  sinusoidally  with  respect  to  the  local  hour  angle.  The  relative  velocity 
of  an  observer  at  an  intermediate  latitude,  moving  from  D  to  E  to  F,  will  also  vary  sinus¬ 
oidally,  but  at  lower  amplitude.  The  amplitude  at  an  intermediate  latitude  is  reduced  by 
the  cosine  of  the  latitude  L  .  This  leads  to  an  expression  of  relative  velocity  of  observer 
with  respect  to  the  moon  of 


V  =  -  <i)  R  cos  L  sin  LHA. 


(3) 


Fig.  3  -  Lines  of  equal  velocity  with  respect  to  moon, 
viewed  from  north  celestial  pole .  Plane  of  paper  is 
earth’s  equatorial  plane.  Parallel  lines  are  inter¬ 
section  of  surface  of  earth  and  planes  normal  to  paper  . 
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Figure  4  illustrates  a  family  of  curves  of  relative  velocity  as  a  function  of  local  hour 
angle.  Figure  5  shows  a  family  of  curves  of  rate  of  change  of  relative  velocity  with  respect 
to  local  hour  angle.  Since  the  Doppler  frequency  shift  is  proportional  to  relative  velocity, 
these  curves  pertain  equally  to  Doppler  shift,  and,  by  substitution  of  velocity  from  Eq.  (3) 
into  the  Doppler  shift  expression,  Eq.  (2), 


Af  =  - 


2  ft 


CO  R  cos  L  sin  LHA. 


Differentiation  with  respect  to  local  hour  angle  yields 


(4) 


d  Af 
d  LHA 


2  f 

- -  CO  R  cos  L  cos  LHA. 

c 


(5) 


Since  c.)  is  equal  to  d  LHA/dt ,  the  time  rate  of  change  of  Af  may  be  found  by  multiplying 
Eq.  (5)  by  If,  at  the  same  time,  Eq.  (4)  is  also  multiplied  by  co,  the  following  pair  of 
parametric  equations  related  by  the  local  hour  results: 


d  Af 
dt 


Af 


2  co^  R  cos  L  cos  LHA 
c 


(6) 


2  f..  R  cos  L  sin  LHA 

CO  Af  =  -  - ^ - (7) 

c 


These  are  the  equations  of  a  sphere  of  radius,  2  f^  co^  R/c.  Dividing  Eq.  (7)  by  Eq.  (6), 
one  obtains 


- ; -  —  wall  • 

Af 

West  longitude  is  determined  by  subtracting  the  local  hour  angle  from  the  Greenwich  hour 
angle;  east  longitude  is  obtained  by  subtracting  the  Greenwich  hour  angle  from  the  local 
hour  angle.  In  symbolic  form. 


Long  W  =  GHA  -  LHA,  and 
Long  E  =  LHA  -  GHA. 


It  may  be  necessary  to  add  360  degrees  to  either  quantity.  The  Greenwich  hour  angle  of 
an  object  is  found  by  subtracting  the  right  ascension  (RA)  of  the  object  from  the  Greenwich 
sidereal  time  (GST)  of  the  observation.  The  Greenwich  hour  angle  of  the  moon  may  be 
found  in  the  Air  Almanac  (or  equivalent  tables). 


When  the  value  of  LHA  has  been  determined  as  above,  Eq.  (4)  may  be  used  to  obtain 
the  latitude  of  the  observer. 


L  =  arccos 


c  Af 

2  f  j  0)  R  sin  LHA 


A  nomograph,  furnishing  a  graphical  procedure  for  the  solution  of  Eqs.  (6)  and(7)  is 
shown  in  Fig.  6.  This,  like  Fig.  3,  is  a  view  of  the  earth  from  the  extended  polar  axis. 
There  is  a  rectangular  reference  system  having  an  ordinate  of  -coAf  and  an  abscissa 
of  -Af .  An  observer  on  the  earth,  capable  of  determining  the  Doppler  frequency  shift  Af , 
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Fig .  4  -  Relative  velocity  with 
respect  to  moon  as  a  function  of 
local  hour  angle 


Fig.  5  -  Relative  acceleration 
with  respect  to  moon  as  a  func¬ 
tion  of  local  hour  angle 


LOCAL  HOUR  ANGLE  (DEGREES) 


LATITUDE  ■ 
0  DEGREES 
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MOONSET 

90 


MOONRISE 

Fig.  6  -  Moon  Doppler  navigation  nomograph 


the  time  rate  of  change  of  Doppler  shift  Af ,  and  the  location  in  the  frequency  domain  of 
the  returned  signal  with  respect  to  the  transmitted  signal,  would  enter  the  nomograph  and 
determine  his  latitude  and  local  hour  angle  with  respect  to  the  moon.  The  nomograph  is 
used  in  the  following  manner: 

1.  Determine  the  maximum  radius  of  the  nomograph  as  r  =  2  ai^  R/c. 

2.  Normalize  the  product  of  earth’s  angular  velocity  and  the  measured  value  of 
Doppler  shift  ojAf/r.  This  equals  the  product  of  cos  L  sin  LHA  and  is  the 
ordinate  of  the  earth  site  on  the  nomograph . 

3.  Normalize  the  measured  value  of  time  rate  of  change  of  Doppler  shift  A’f/r  . 

This  equals  the  product  of  cos  L  cos  LHA  and  is  the  abscissa  of  the  earth 
site  on  the  nomograph. 

A  point  having  a  given  ordinate  and  abscissa  is  shown  in  Fig.  6. 

This  particular  site  has  a  latitude  of  30  degrees  N  and  a  moon  local  hour  angle  of 
25  degrees.  The  radial  lines  from  the  origin  in  Fig.  6  represent  the  local  hour  angle  of 
the  moon  from  the  site.  The  fourth  quadrant,  in  the  mathematical  sense,  represents  local 
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hour  angles  from  moonrise  (LHA  =  270  degrees)  to  the  local  meridian  (lha  =  0  degree). 
The  first  quadrant  represents  local  hour  angles  from  the  local  meridian  to  moonset 
(LHA  =  90  degrees).  The  second  and  third  quadrants  represent  local  hour  angles  from 
90  to  270  degrees  and  are  required  when  the  declination  of  the  moon  varies  from  zero  degre 

The  concentric  circles  are  parallels  of  latitude  for  the  northern  hemisphere  as  shown 
in  the  nomograph.  A  view  of  the  nomograph  from  below  the  plane  of  the  paper  would 
represent  the  southern  hemisphere  (in  other  words,  there  is  a  hemispheric  ambiguity  in 
the  nomograph  and  in  the  solution  of  the  equations  on  which  it  is  based). 

A  simplified  system  has  been  assumed  for  purposes  of  clarity.  The  existing  situation 
is  more  complex  than  the  one  just  described.  The  moon’s  orbital  velocity,  changing  dis¬ 
tance  from  the  earth,  and  varying  declination  enter  the  calculation  of  relative  velocity. 

In  addition,  the  transmitter  may  be  widely  separated  from  the  receiver  (bistatic  operation). 
Further  discussion  of  these  effects  appears  in  Appendix  D.  The  principles  described 
above  are  not  changed  even  when  the  complexity  of  the  existing  earth-moon  system  is  con¬ 
sidered.  A  system  of  measurements  based  on  these  principles  would  give  a  continuous 
flow  of  position  data  during  the  time  from  moonrise  to  moonset.  Such  a  system  could  be 
used  as  the  base  of  a  terrestrial  navigation  system,  independent  of  weather  (which  affects 
celestial  navigation  by  optical  means). 


PROBABLE  ERROR  IN  POSITION  DETERMINATION 

When  the  value  of  a  property  is  to  be  computed  from  observations  of  other  properties, 
it  is  natural  to  inquire  about  the  exactness  with  which  each  observation  should  be  performed 
to  produce  a  given  degree  of  precision  in  the  result.  Precision  expresses  the  degree  of 
concordance  or  consistency  within  a  set  of  observed  values.  The  best  representation  of  the 
degree  of  consistency  in  a  set  of  discordant  measurements  of  a  quantity  which  were  per¬ 
formed  with  equal  exactness  is  the  probable  error.  The  probable  error  of  a  quantity  may 
be  defined  as  that  error  which  is  just  as  likely  to  be  exceeded  as  not  to  be  exceeded.  The 
probable  error  of  a  computed  result  is  derived  in  Appendix  C. 

The  position  of  a  site  on  the  surface  of  the  earth  is  calculated  from  measurements  of 
Doppler  frequency  shift  and  time  rate  of  change  of  Doppler  shift,  which  in  turn,  are 
dependent  on  the  values  of  relative  velocity  and  acceleration  of  the  site  with  respect  to  the 
moon.  If  the  probable  error  in  position  is  desired,  it  is  necessary  to  determine  the 
probable  error  in  both  the  Doppler  frequency  shift  and  the  time  rate  of  change  of  Doppler 
shift.  This  is  done  in  Appendix  E. 

The  equation  of  relative  velocity  of  a  site  on  the  earth  with  respect  to  the  moon  is 
derived  in  Appendix  D.  The  derivation  of  an  equation  describing  the  relative  acceleration 
of  the  site  with  respect  to  the  moon  from  the  complete  velocity  equation  is  a  formidable 
task,  however,  it  is  illuminating  to  develop  the  probable  errors  in  relative  velocity  and 
acceleration  from  simplified  velocity  and  acceleration  equations  based  on  earth  rotation 
alone.  The  order  of  magnitude  of  errors  in  position  obtained  in  this  case  may  be  said  to 
represent  the  minimum  probable  errors  to  be  expected.  This  analysis  should  indicate  the 
degree  of  usefulness  of  the  Doppler  method  of  navigation. 

The  probable  errors  in  relative  velocity  and  acceleration  due  to  the  rotation  of  the 
earth  are  derived  in  Appendix  E,  and  the  position  errors  in  latitude  and  longitude  are 
derived  in  Appendix  F.  The  errors  in  position  are  functions  of  position  as  well  as  time. 

The  extreme  values  of  position  errors  in  latitude,  as  a  function  of  latitude,  for  local  hour 
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angles  of  0  and  90  degrees  are  shown  in  Fig.  7.  The  longitude  errors,  for  local  hour 
angles  of  0  and  90  degrees  as  a  function  of  latitude  are  shown  in  Fig.  8.  The  “assumed 
geodetic  error  in  position”  is  based  on  the  assumption  that  a  position  on  the  surface  of  the 
earth  may  be  determined  to  within  plus  or  minus  ten  feet  in  three  mutually  perpendicular 
directions  by  means  of  geodetic  and  gravimetric  surveys. 


Fig.  7  -  Latitude  error  as  a  function  of  latitude 


A  plot  of  contours  of  constant  probable  error  in  latitude  (in  kilometers),  when  the 
geographical  point  of  the  moon  is  at  latitude  0  degree  and  longitude  0  degree,  is  shown 
in  Fig.  9.  This  plot  shows  that  navigation  to  within  line  of  sight  of  a  given  position  may 
be  accomplished  over  70  percent  of  the  hemisphere  visible  to  the  moon.  This  a  clear 
indication  of  the  degree  of  usefulness  of  moon  Doppler  navigation. 


TRANSMITTER  FREQUENCY  STABILITY 


If  the  Doppler  frequency  shift  is  determined  by  beating  the  transmitter  frequency  with 
the  echo,  the  quantity  actually  measured  at  the  earth  site  is 


Af' 


(' 


+ 


- 
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Fig.  8  -  Longitude  error  as  a  function 
of  latitude  (declination  =  0  degrees) 


where 

f  ^  is  the  transmitter  frequency  at  the  time  of  transmission 
f  ^ '  is  the  transmitter  frequency  at  the  time  of  reception  of  the  echo 
V  is  the  relative  velocity  of  site  and  moon,  and 
c  is  the  velocity  of  propagation. 

The  measured  value  of  Doppler  shift  is  a  function  of  four  variables,  Af '  =  G(f^,  v,c,  f^'), 
and  the  probable  error  in  the  measured  value  of  Doppler  shift  is  given  by 


^Af 


or 


■  ■  /(ft'  •(‘a"  ••)'  •  te'  -)■  •  («7  ■")■ 


«Af ' 


ec^  +  Eft, 


2  . 


Fig.  9  -  Contours  of  constant  probable  error  in  latitude  (in  kilometers)  when  the  geographical  point  of  the  moon  is 

at  latitude  0  degrees  and  longitude  0  degrees 
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Since  v/c  «  1,  this  expression  reduces  to 


^Af 


,  /  'vV 

*v— ; 


An  examination  of  the  range  of  values  taken  on  by  the  probable  error  in  relative  velocity 
(Appendix  E)  gives  these  values  as  extremes  at  low  and  high  latitudes  respectively: 

e^(max)  .=  6.4  x  10-^  meters/sec 

e^(min)  .=  2.2  x  10"^  meters/sec. 

It  can  be  seen  from  the  expression  for  the  probable  error  in  measured  Doppler  shift  e^j/, 
that  the  ratio  of  probable  error  in  transmitter  frequency  to  transmitter  frequency  should 
be  an  order  less  than  the  ratio  of  probable  error  in  relative  velocity  to  the  velocity  of 
propagation,  if  the  probable  error  in  measured  Doppler  shift  is  to  be  primarily  dependent 
on  the  latter  ratio.  If  this  is  the  case,  the  requirement  of  transmitter  frequency  stability 
in  this  system  may  be  expressed  as 

1  «ft  1 

-  <  - — —  <  -  , 

10^3  "  ft 


however,  at  the  higher  latitudes  (where  the  frequency  stability  requirements  approach 
1  part  in  10*^)  the  error  in  position  diminishes,  and  for  a  specified  error  in  position  the 
stability  requirements  may  be  relaxed. 


LIBRATION  EFFECT 

If  the  surface  of  the  moon  is  rough  and  scatters  radiation  in  all  directions,  there  will 
be  an  increase  in  the  directivity  of  the  scattered  radiation  and  of  the  echo  cross  section 
compared  to  that  from  a  smooth  moon.  Energy  may  then  be  returned  from  the  whole  visible 
disk  and  not  only  from  the  first  Fresnel  zone.  This  introduces,  if  the  moon  is  rotating, 
a  Doppler  spread  (equal  to  4  vA  )  in  the  echo  spectrum  corresponding  to  plus  or  minus 
the  velocity  of  approach  v  of  the  limb  relative  to  the  center  (2).  The  angular  velocity  of 
libration  of  the  moon  reaches  a  maximum  of  about  3  degrees  per  day,  hence,  the  velocity 
of  approach  of  the  limb  is  of  the  order  of  one  meter  per  second.  Although  the  position 
errors  in  latitude  and  longitude  are  shown  (Appendix  F)  to  be  independent  of  transmitter 
frequency  when  the  Doppler  shift  due  to  earth  rotation  alone  is  employed,  the  libration 
effect,  as  a  source  of  error,  would  increase  with  increasing  frequency.  However,  good 
evidence  has  been  obtained  recently  which  indicates  a  reflecting  zone  much  smaller  than 
the  visible  disk  of  the  moon.  This  would  greatly  reduce  the  width  of  the  echo  spectrum. 


EXPERIMENT 

Two  experiments  were  performed,  one  on  May  13,  1956,  the  other  on  July  7,  1956,  at 
field  sites  of  the  laboratory  for  the  purpose  of  measuring  the  Doppler  frequency  shift 
undergone  by  a  301-Mc  transmission  from  earth  to  moon  to  earth.  Each  experiment  had 
a  duration  of  approximately  thirty  minutes  because  of  the  limited  aspect  of  the  transmitting 
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antenna.  The  transmitting  site  comprised  a  10-kw,  cw,  301 -Me  transmitter  and  a  fixed 
parabolic  antenna  having  a  gain  of  37  db  at  this  frequency.  The  physical  aperture  of  the 
antenna  had  dimensions  of  220  ft  by  263  ft  and  is  shown  in  Fig.  10.  Limited  steering  was 
accomplished  by  moving  the  antenna-feed  system  on  either  side  of  the  focus  of  the  parabo¬ 
loid  within  a  region  of  low  aberration.  The  receiving  site  included  a  receiver  having  a 
noise  figure  of  6  db  fed  by  a  25-ft  paraboloid  with  a  gain  of  25  db.  The  receiving  antenna 
is  shown  in  Fig.  11.  The  system  parameters  were: 


Transmitter  power  10  kw 

+40  dbw 

Transmitting  antenna  gain 

+  37  db 

Receiving  antenna  gain 

+25  db 

Path  attenuation 

-258  db 

Available  signal  power  at  receiver 

-156  dbw 

-156  dbw 

Noise  power/cycle  (kT)  4. 14  X10“^^  watts 

-204  dbw 

1-kc  bandwidth  10®  cps 

+  30  db 

Noise  figure  of  receiver 

+  6  db 

Noise  power/kc  bandwidth  (signal  power 
for  S/N  ratio  of  Odb) 

-168  dbw 

-168  dbw 

Estimated  signal  to  noise  ratio  -156  -  (-168) 

+12  db 
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Fig  .11-  Receiving  antenna 


The  Doppler  frequency  shift  was  obtained  by  beating  the  returning  signal  with  the  ground 
wave  arriving  at  the  receiving  site  from  the  transmitting  site  25  miles  distant,  Fig.  12. 

An  audio  oscillator  was  tuned  to  the  beat  frequency  by  visual  observation  of  Lissajous 
figures  on  an  oscilloscope.  An  electronic  counter  was  employed  to  measure  the  audio 
oscillator  frequency.  The  measured  values  of  the  Doppler  shift,  as  a  function  of  time, 
are  compared  with  the  computed-value  curves  in  Figs.  13  and  14. 

The  absolute  value  of  the  average  error  (or  deviation  of  measured  values  from  the 
computed  curve)  is  approximately  4  cps.  The  computed  value  of  probable  error  in  Doppler 
shift  (under  the  experimental  conditions)  is  0. 1  cps.  The  application  of  “the  method  of  least 
squares”  to  find  the  best  linear  representation  of  the  experimental  values,  results  in  a 
curve  which  agrees  with  the  computed  curve  within  plus  or  minus  0. 1  cps. 


CONCLUSIONS 

The  location  of  a  fixed  site  on  the  earth  may  be  determined  within  plus  or  minus  10  feet 
in  three  mutually  perpendicular  directions  by  means  of  geodetic  and  gravimetric  surveys. 
These  surveys  are  performed  with  great  care  and  require  considerable  data.  It  is  logical 
to  expect  the  precision  of  a  navigational  system  to  be  proportional  to  the  quantity  and 
quality  of  data  employed.  It  has  been  the  purpose  of  this  report  to  examine  and  evaluate 
the  quality  of  available  data  employed  in  a  moon  Doppler  navigation  system. 

The  fundamental  quantities  employed  to  determine  position  on  the  surface  of  the  earth, 
in  this  system,  are  the  Doppler  frequency  shift  and  the  time  rate  of  change  of  Doppler  shift 
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Fig.  13  -  Experimentadly  determined 
Doppler  frequency  shift  as  a  function 
of  time 
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Fig.  14  -  Experimentally  determined  Doppler 
frequency  shift  as  a  function  of  time 


undergone  by  a  radio -frequency  signal  after  it  has  been  reflected  from  the  moon  and  received 
on  earth.  The  following  conclusions  have  been  deduced  from  the  study  presented  in  this  report: 

1.  There  is  essentially  a  north-south  hemispheric  ambiguity  in  Doppler 
position  determination. 

2.  Latitude  and  longitude  position  errors  inherent  in  a  moon  Doppler  navigation 
system  are  functions  of  position  as  well  as  time. 

3.  Position  errors  increase  with  decreasing  latitude  (latitude  position  errors 
range  from  a  few  irieters  at  90-degree  latitude  to  100  kilometers  at  0-degree 
latitude;  and  longitude  position  errors  range  from  a  few  meters  at  90 -degree 
latitude  to  about  2  kilometers  at  0-degree  latitude),  however,  navigation  to 
within  line  of  sight  of  a  given  position  (10  kilometers)  may  be  accomplished 
over  70  percent  of  the  hemisphere  illuminated  by  the  moon  at  a  given  time. 

4.  Position  errors  are  independent  of  the  transmitter  frequency  employed. 

5.  The  probable  error  in  transmitter  frequency  must  be  within  the  order  of 
10^*  if  the  full  capabilities  of  the  system  are  to  be  realized. 
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This  study  is  based  on  an  idealized  earth-moon  system  where  the  relative  motion 
between  earth  site  and  the  moon  is  that  due  to  rotation  of  the  earth  alone.  However,  an 
indication  that  the  position  errors  described  are  of  the  right  order  of  magnitude  is  the 
fact  that  the  application  of  the  method  of  least  squares  to  the  limited  experimental  Doppler 
measurements  yields  results  which  agree  with  the  calculated  values  of  Doppler  frequency 
shift  within  the  computed  value  of  the  probable  error. 

It  is  apparent  that  the  errors  in  position  inherent  in  the  Doppler  method  of  position 
determination  are  large  at  low  latitudes,  however,  the  errors  diminish  at  high  latitudes, 
where  some  navigational  systems  are  relatively  inadequate.  Furthermore,  an  angle¬ 
measuring  method  (radio  sextant)  might  be  employed  to  complement  the  Doppler  method 
at  low  latitudes. 
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APPENDIX  A 

List  of  Magnitudes  and  Probable  Errors  of  Parameters 


c  =  2.997893  X  10®  meters/second* 
e^,  =  5  X  10®  meters/second 
R  =  6.378260  X  10®  meterst 
Ejj  =  1  X  10®  meters 
p  -  1.748  X  10"^  radians 
p  =  2.042  X  10"^  radians 
=  5  X  10"®  radians 
dp/dt  =  -5  X  10"^®  radians/second 
=  8.1  X10"‘®  radians/second 
d2p/dt2  =  -7.  lx  10'^®  radians/second  ^ 

=  1.3X10"‘*  radians/second® 
sin  p  =  0.017813 
cos  p  =  0.99984 

=  1  X  10“®  seconds 
5  =  20°  23' 


(velocity  of  propagation) 

(probable  error  in  c ) 

(equatorial  radius  of  earth) 

(probable  error  in  R) 

(horizontal  parallax  angle  of  moon,  May  13) 
(horizontal  parallax,  July  7) 

(probable  error  inp) 

(May  13,  1956) 

(July  7,  1956) 

(May  13,  1956) 

(July  7,  1956) 

(probable  error  in  time) 

(declination  of  moon.  May  13,  1956  at  2230  UT) 
(declination  of  moon,  July  7,  1956  at  1900  UT) 
(probable  error  in  8) 

(May  13,  1956) 

(July  7,  1956) 


*Bol  and  Hansen,  “Microwave  Cavity  Determination,  ”  Stanford  University,  1950 
tArmy  Map  Service,  1954 
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L  =  38. 545  degrees 

(latitude  of  transmitter) 

6l  =  5  X  10"^  radians 

(probable  error  in  site  latitude,  if 
surveyed  by  geodetic  and  gravimetric  means) 

cos  L  =  0.78212 

sin  L  =  0. 62303 

cj  =  7  X  10"*  radians/second 

(angular  velocity  of  earth  with  respect  to 
moon) 

e  =  1  X  10”® 

eo 

(probable  error  in  oj ) 

LHA  =  0.557  radians 

(local  hour  angle  of  moon) 

^LHA  =  1  ^  10”*  radians 

(probable  error  in  LHA) 

sin  LHA  =  0.55194 

cos  LHA  0.83389 

♦  *  * 


APPENDIX  B 

Derivation  of  Doppler  Frequency  Shift  and 
Time  Rate  of  Change  of  Doppler  Shift 


If  a  receiver  is  in  motion  toward  a  stationary  transmitter,  the  received  frequency 
will  be 

f,  =  f,  (1  +  v/c)  (Bl) 

and  if  a  transmitter  is  in  motion  toward  a  stationary  receiver,  the  received  frequency 
will  be 


where 

is  the  frequency  of  the  received  signal 
is  the  frequency  of  the  transmitted  signal 

c  is  the  velocity  of  propagation  in  the  medium  between 
transmitter  and  receiver 

V  is  the  relative  velocity  of  transmitter  and  receiver. 

If  the  moon  is  in  motion  relative  to  a  site  on  the  earth,  and  a  signal  is  sent  from  the 

site  to  the  moon  and  there  reflected,  the  frequency  of  the  signal  received  at  the  earth  site 

will  be  1 

^  =  ft  (1+''/c)(yt^)-  (B3) 

If  the  quantity,  (1  -v/c)-i  ,  is  expressed  as  a  series  and  multiplied  by  (1  +v/c),  the  result  is 

fj,  =  (1  +  2v/c  +  2v2/c^  +  2v^/c3  +  ...)  (B4) 

which,  if  terms  of  higher  order  than  the  first  power  are  neglected,  gives 

fj.  =  fj,  (l+2v/c),  v/c  <  1.  (B5) 

The  Doppler  frequency  shift  is  defined  as  the  difference  between  the  transmitted  and 
received  frequencies. 

Af  =  f,.  -  X  (2v/c).  (B6) 

The  time  rate  of  change  of  Doppler  frequency  shift  is  simply 

d  Af  _  2  ft  dv 
dt  c  dt 

*  *  * 
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APPENDIX  C 

Probable  Error  of  a  Computed  Result 


If  y  =  f(x)  ,  an  error  in  x  will  produce  an  error  in  y .  If  the  errors  in  x  and  y  are 
denoted  by  Ax  and  a  y,  the  ratio  of  errors  is  Ay/Ax,  and  as  Ax  approaches  zero  Ay/Ax 
approaches  dy/dx  and 


dy  = 


(Cl) 


The  differentials  dy  and  dx  represent  errors  in  y  and  x  respectively.  When 
y  =  f(x,z,v, . . . ), 


Squaring  both  sides  of  Eq.  (C2)  gives 


(C2) 


(C3) 


If  the  errors  are  symmetrically  distributed  about  the  mean,  the  cross  products  will 
cancel  and  only  the  squared  terms  will  remain  provided  the  number  of  sets  is  sufficiently 
large.  Therefore, 
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APPENuia  D 

Derivation  of  Relative  Velocity  Between  Earth  Site  and  Moon  Center 
A  view  of  the  earth -moon  system  is  shown  in  Fig.  D1 


Fig.  D1  -  The  earth-moon  system 


where 

A  is  a  terrestrial  site  at  latitude  L 
0  is  the  earth  center 
X  is  the  moon  center 
S  is  the  declination  of  the  moon 
R  is  the  radius  of  the  earth 

HA  is  the  distance  from  the  site  to  the  earth’s  polar 
axis  and  equals  R  cos  L 

OH  is  the  distance  from  the  earth’s  equatorial  plane 
to  the  plane  of  latitude  L  and  equals  R  sin  L 

LHA  is  the  local  hour  angle  of  the  moon  measured  from  site  A 

OX  is  the  distance  from  earth  center  to  moon  center 

CB  is  the  distance  from  the  earth’s  equatorial  plane  to  the 
plane  in  which  the  moon  appears  to  move  due  to  the 
rotation  of  the  earth  and  equals  OX  sin  s 

AX  is  the  distance  from  the  site  to  moon  center 

HL  is  the  distance  from  the  earth’s  polar  axis  to  the  moon  center 
and  equals  OX  cos  s. 
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The  equations  relating  the  distance  AX  to  the  other  quantities  are  constructed  as 
follows: 

(AX)2  =  (AL)2  +  (HB)2  (D1) 

(AL)2  =  (HA)2  +  (HL)2  -  2  (HA)  (HL)  cos  LHA 

=  r2  cos2  L  +  (0X)2  cos2  8  -  2R  (OX)  cos  L  cos  8  cos  LHA  (D2) 

HB  =  OB  -  OH  =  OX  sin  8  -  R  sin  L.  •  (D3) 

The  substitution  of  Eqs.  (D2)and(D3)  into  (Dl)  gives 

(AX)^  ,=  r2  cos2  L  +  (OX)2  cos2  8  -  2R  (OX)  cos  L  cos  8  cos  LHA 
+  (0X)2  sin2  8  +  R2  sin2  L  -  2R  (OX)  sin  8  sin  L. 

This  simplifies  to 

(AX)2  =  r2  +  (0X)2  -  2R  (OX)  (cosLcos  8cos  LHA  +  sinLsinS).  (D4) 

The  distance  AX  is  a  function  of  ox  and  angles,  LHA  and  s,  hence,  the  relative 
velocity  between  A  and  X  is 

..  _  d  (AX)  _  r3  (AX)n  fd  LHA-j  (AX)!  fd  (OX)!  ,  fB  (AX)1  fdSl 

The  individual  terms  of  Eq.  (D5)  are  the  following: 


B  (AX) 
B  LHA 


_ R  (OX)  (cos  L  cos  8  sin  LHA) _ 

^r2  +  ( OX ) 2  -  2R  (OX)  (cos  Lcos  8  cos  LHA+  sinLsinS) 
d  LHA 


(D6) 

(D7) 


The  angular  velocity  ai  is  equal  to  the  difference  between  the  earth’s  angular  velocity 
(sidereal)  and  the  angular  velocity  of  the  moon  in  its  orbit,  or  in  other  words,  w  is  the 
earth’s  angular  velocity  with  respect  to  a  line  connecting  earth  and  moon  centers. 


B  (AX) 
B  (OX) 


OX  -  R(cosLcos  8  cos  LHA+  sin  L  sin  8) 

^R2  +  (OX)  2  -  2R  (OX)  (cosLcos  8  cos  LHA  +  sinL  sin  8) 


(D8) 


(D9) 

tDlO) 
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where  p  is  the  horizontal  parallax  angle  of  the  moon.  This  angle  may  be  defined  as  the 
angle  subtended  by  the  earth’s  radius  observed  from  the  center  of  the  moon. 


B  (AX) 
B  8 


y  sin^ 


-R  (sin  L  cos  8  -cos  L  cos  LHA  sin  S)  COX') _ _ 

p  +  1  -  2  sin  p  (cos  L  cos  S  cos  LHA  + sin  L  sin  S) 


(Dll) 


The  substitution  of  Eqs.  (D6),  (D7),  (D8),  (D9),  and^Dll)  into  (D5)  gives  the  following 
expression  for  relative  velocity  between  the  site  A  and  the  moon’s  center  X 


R  cos  L  cos  8  sin  LHA 


(cos  L  cos  8  cos  LHA  +  sin  L  sin  8)  R  cot  p 


d  t 


R  cot  p  CSC  p 


dt 


D 


(sin  L  COS  8  -  cos  L  cos  LHA  sin  8)  R 


d8 

dt 


D 


(D12) 


where  D  is  ysin^  p  -  2p  (cos  L  cos  8  cos  LHA  +  sin  L  sin  8)  +  1  .  The  quantity  sin  p 
may  be  replaced  by  p  ,  since  p  is  of  the  order  of  0.  02  radian,  and  the  absolute  value  of 
the  error  in  making  the  substitution  is  less  than  p^/d! ,  which  is  less  than  6  parts  in  a 
million.  Similarly,  cot  p  may  be  replaced  by  1/p  with  an  error  of  the  same  order  of 
magnitude.  Relative  velocity  is  then  expressed  by 

V  =  —  ^  —  I  o)  cos  L  cos  8  sin  LHA 

y p'^  -  2p  (cos  L  cos  8  cos  LHA  +  sinL  sin  8)  +  1  L 


^co 


+  ( cos  L  cos  8  cos  LHA  +  sin  L  sin  8 


-  (sin  L  cos  8  -  cos  L  cos  LHA  sin  8) 


PJ  P 

ill 

dt  J’ 


dp 

dt 


(D13) 


The  term  involving  the  time  rate  of  change  of  declination  is  significant  only  when  declina¬ 
tion  approaches  zero,  therefore,  the  relative  velocity  equation  may  be  simplified  to  be 


y p^  -  2p  (cos  L  cos  8  cos  LHA+  sin  L  sin  8) 


c 

1  - 


cos  L  cos  8  sin  LHA 


+  (  cos  L  cos  8  cos  LHA  +  sin  L  sin  8 


p/  p 


p/  p  dt 


-  sin  L 


d8 

dt 


(D14) 
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The  first  (rotational)  term  represents  a  sinusoidal  variation  having  a  peak  value  of 
approximately  448  meters  per  second  and  a  period  of  24  hours,  at  a  latitude  of  zero  degree. 
Superimposed  on  this  variation  is  that  due  to  the  change  in  orbital  distance  (second  term), 
which  has  a  peak  value  of  about  51  meters  per  second  and  a  period  of  about  28  days.  The 
third  term  (due  to  changes  in  declination)  has  a  peak  value  of  6  meters  per  second  and  a 
period  of  about  28  days.  The  maximum  relative  contributions  of  the  orbital  and  declina¬ 
tion  terms  are  about  11  percent  and  1  percent,  respectively,  of  that  due  to  rotation.  It  is 
not  valid  to  compute  probable  errors  in  relative  velocity  in  the  accustomed  manner  from 
this  equation,  since  an  error  in  velocity  is  not  only  a  function  of  position  but  of  time. 
Although  the  relative  contribution  of  declination  changes  to  velocity  is  but  1  percent  of  that 
due  to  rotation,  there  are  times  when  the  declination  change  is  the  sole  source  of  relative 
velocity  between  earth  site  and  moon.  The  same  situation  applies  to  the  contribution  of 
orbital  distance  changes  to  relative  velocity. 

Another  point  of  interest  is  the  denominator  of  the  coefficient,  or  multiplying  factor, 
of  the  equation.  The  extreme  values  of  the  denominator  are  1.0002,  and  0.98.  If  the 
denominator  is  assumed  to  be  unity,  as  is  sometimes  done,  an  error  of  2  percent  may  be 
realized  in  the  computed  value  of  relative  velocity. 

The  derivation  of  an  equation  representing  the  acceleration  of  an  earth  site  relative 
to  the  moon  from  Eq.  (D14)  is  a  formidable  task,  however,  for  the  purpose  of  examining 
an  extremely  idealized  case,  the  velocity  equation  will  be  simplified  to  be  v  =  ojRcos  L  cos  s 
sin  LHA,  The  equation  for  acceleration  is  then  simply 

dv 

-^  =  V  =  or  R  cos  L  cos  8  cos  LHA.  (D15) 

The  value  of  acceleration  given  by  this  equation  is  due  to  the  rotation  of  the  earth  alone, 
but  values  of  probable  errors  derived  from  this  and  the  simplified  velocity  equation  will 
be  enlightening  and  give  an  idea  of  the  order  of  m^nitude  of  errors  in  position  to  be 
expected  in  a  navigation  system  based  on  the  measurement  of  relative  velocity  and  accelera¬ 
tion.  This  is  done  in  Appendixes  E  and  F. 


*  *  * 
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APPENDIX  E 

Derivation  of  Probable  Error  in  Computed  Values  of  Doppler 
Frequency  Shift  and  Time  Rate  of  Change  of  Doppler  Shift 


The  Doppler  frequency  shift  (Appendix  B)  is  given  by  Af  =  2  f^v/c,  therefore,  the 
probable  error  in  Af  (Appendix  C)  is 


fBAfV  ^  p  Afl  ^  r^AfT 

.3  f  J  "ft  [3  V  J  ^  [ac. 


Performing  the  indicated  differentiation  and  simplifying  gives 


The  time  rate  of  change  of  Doppler  shift  is  given  by  Af  =  2  f^/c  v  ,  and  an  analogous 
operation  gives 


€aV  -  Af 


It  is  of  interest  to  note  that  the  quantities  actually  measured  by  an  observer  at  the 
earth  site  are  not  Af  and  A'f.  The  quantities  measured  are,  if  Af  is  determined  by  beat¬ 
ing  the  transmitter  frequency  with  the  echo. 


Af'  =  f,  -  ft(t  +  At) 


The  quantity  f  ^  (Appendix  B)  is 


fr  =  ft  + 


2  f^  V 


hence 


Af'  =  Af  -  Af» 


where  Af^  is  the  change  in  transmitter  frequency  during  the  transit  time  of  the  signal 
from  earth  to  moon  to  earth. 


The  relative  velocity  of  a  site  on  the  earth  with  respect  to  the  moon  due  to  rotation 
of  the  earth  (Appendix  D)  has  been  shown  to  be  v  =  coR  cos  L  cos  5  sin  LHA  .  The  probable 
error  in  relative  velocity  is 


(S]^ 


[51]'  ,  [51]'  ,  : 

LBtJ  ‘  LSojJ 
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The  partial  derivatives  are: 


3v 

-  =  o)  cos  L  cos  S  sin  LHA 

3R 


R  cos  L  cos  S  cos  LHA 


—  =  R  cos  L  cos  S  (sin  LHA  + LHA  cos  LHA) 

Bc<j 


3v  .  ,  . 

—  =  -cijR  sin  S  cos  L  sin  LHA 
BS 


~  =  -toRcos  S  sin  L  sin  LHA*  (ElO) 

BL 

The  relative  acceleration  of  a  site  on  earth  with  respect  to  the  moon  due  to  rotation 
of  the  earth  (Appendix  D)  is  v  =  R  cos  L  cos  8  cos  LHA.  The  probable  error  in 
acceleration  is 

-(er  --  [£r 

*  ,.n, 


The  partial  derivatives  are: 


cos  L  cos  5  cos  LHA 
BR 


(E12) 


—  =  -  R  cos  L  cos  8  sin  LHA 


(E13) 


— —  =  R  cos  L  cos  S  (2ajcos  LHA  -  sin  LHA) 
Bo; 


{E14) 


3  V  , 

—  =  -  0}^  R  cos  L  sin  8  cos  LHA 
BS 


{E15) 


-  co^  R  sin  L  cos  8  cos  LHA. 


(E16) 
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APPENDIX  F 

Derivation  of  Probable  Errors  in  Position 


The  general  equation  of  a  sphere  of  radius  r  in  a  rectangular  coordinate  system  is 


x2  +  y2  +  z2  .=  r2 


where 


X  =  r  cos  cos  6 
y  =  r  cos  sin  6 
z  =  r  sin  ^ 


^  . '  y 

6  -  arctan — 

X 


(f>  =  arccos 


i/^ 


The  sphere  described  by  these  equations  is  shown  in  Fig.  FI. 


Fig.  FI 
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The  error  in  position  on  the  surface  of  the  sphere  is  Tt4>  the  4>  direction,  and 
r  cos  4>  A0  in  the  e  direction.  The  angles  4>  and  e  are  functions  of  x  and  y ,  therefore, 
the  incremental  angles  and  AS  are  found  in  the  following  manner; 


cos  = 


cos  ((^  +  A</)) 


^(x  +Ax)2  +  (y+Ay)^ 


/x2  +  y2  2  (x  Ax  +  y  Ay)  Ax^  +  Ay  ^ 

-  +  -  +  - . 

r2  r2  r2 

If  higher  orders  than  the  first  of  Ax  and  Ay  are  neglected  and  substitutions  for  x  and  y 
are  made,  the  result  is 

r  T  2  cos  '  , 

cos  cfi  cos  A<^  “  sin  (j>  sin  ^4>  -  t/  cos-^  (p  +  -  (cos  0  Ax  +  sin  0  A  y  ) 


sin  <p  sin  A^  =  [cos  <p]  cos  -  1  / 1  + 


2  (cos  gAx  +  sin  0  Ay) 
r  cos  <p 


When  A<^>  approaches  zero,  Eq.  (FI)  becomes 


sin  <p  dp  -  [cos  1  -  1 /l  + 


2  (cos  0  dx  +  sin  0  dy) 
r  cos  p 


and  since  Yl  +  A  %  l  +  A/2  when  A  is  very  small. 


sin  p  dp  =  [cos  p] 


-  (cos  0  I 

T 


dx  +  sin  S  dy  ) 
cos  p 


-  (cos  0dx  +  sin  0  dy) 

dp  -  - : -  •  (F. 

r  sin  p 

This  is  the  angular  error  in  the  p  direction  as  a  function  of  x ,  y,  and  9 ,  and  may  be 
used  when  p  >  hp .  However,  when  p  becomes  small  or  approaches  zero,  Eq.  (FI) 
becomes 


cos  hp  =  1/  1  + 


2  (cos  0  Ax  +  sin0  Ay) 


cos2  A<^  =  1  + —  (cos  0  Ax  +  sing  Ay) 
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sin2  A(^  =  —  (cos  (9  Ax  +  sin  (9  Ay) 


sin  A(^  =  Ai^ 


(cos  0  Ax  +  sin  6  Ay) 


This  is  the  expression  for  the  determination  of  the  angular  error  in  4>  for  values  of 


near  zero. 


The  angular  error  in  the  B  direction,  AB,  is  determined  as  follows: 


tan  B  =  — 

X 


tan  (B+AB)  - 


y  +  Ay 
X  +  Ax 


tan  B  +  tan  AB  _  y  +  Ay 
1  -  tan  B  tan  AB  x  +  Ax 


and  replacing  tan  AB  by  aB, 


which  reduces  to 


(A(9+tan  B)  (x+Ax)  =  (1-A0  tan  6)  (y  +  Ay) 

_  y  +  Ay  -  (x  +Ax)  tan  B 
X  +  Ax  +  (y  +Ay)  tan  6 

r  cos  <i  sin  0  +  Ay  -  ~  (r  cos  4>  cos  9  + Ax) 

^  cos  B _ 

r  cos  (b  cos  B  +  Ax  +  T  (r  cos  <b  sin  B  +  Ay ) 

^  H 


cos  B  Ay  -  sin  B  Ax _ 

r  cos  b  +  cos  B  Ax  +  sin  B  Ay 


In  the  general  case,  where  4>  does  not  ^proach  90  degrees,  AB  may  be  represented  by 

cos  B  Ay  -  sin  B  Ax  /--s 

AB  -  - - 

r  cos  (f> 

The  limiting  value  of  AB  when  4>  approaches  90  degrees  is 

lim  .n  =  cos  B  Ay  -  sin  B  Ax_  ^ 
cos  6>  Ax  +  Sin  B  Ay 

The  equations  describing  the  relative  velocity  and  acceleration  of  an  earth  site  with 
respect  to  the  moon  due  only  to  the  rotation  of  the  earth  (Appendix  D)  are 

V  =  o)  R  cos  S  cos  L  sin  LHA 


V  =  R  cos  8  cos  L  cos  LHA . 
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The  Doppler  frequency  shift  and  the  time  rate  of  change  of  Doppler  shift  (Appendix  B)  are 
given  by 


and 


When  Af 


Af 


Af 


is  multiplied  by 


2  fj,  V  2  fj  coR  cos  8  cos  L  sin  LHA 
c  c 

2  fj.  V  2  fj.  a)2R  cos  8  cos  L  cos  LHA 
c  c 

oi,  the  equations  for&j  Af  and  Af  describe  a  sphere 


where 

X  =  Af 
y  =  6j  Af 

Ax  =  e  •  (probable  error  in  the  time  rate  of  change  of  Doppler  shift) 

Af 

Ay  =  e  =  i/ 6.  2  +  /^f2  g  2  (probable  error  in  the  product  of  co  and  Af) 

o)  Ijit  r  Ar 

4>  =  L  (latitude) 

0  =  LHA  (local  hour  angle) 

2  f*  R  cos  8 

r  = - £ -  . 

c 


The  substitution  of  these  quantities  in  Eqs.  (F3),  (F4),  (F5),  and  (F7)  gives  expressions 
for  the  angular  errors  in  latitude  and  longitude. 


The  probable  position  error  in  latitude  (in  meters)  is  given  by  the  product  of  the 
angular  error  in  latitude  and  the  radius  of  the  earth.  Using  Eq.  (F3) 


R  dL 


cos  LHA  +  sin  LHA 


+  Af2 


2  fj,  cos  8  sin  L 


(FS) 


which  approaches  Eq.  (F4)  as  a  limit,  when  latitude  approaches  zero, 


R  AL  = 


-  c  R  cos  LHA  +  sin  LHA  y/  o?'  +  Af^ 


(F9) 


fj.  ct)^  cos  8 


The  probable  position  error  in  longitude  (in  meters)  is  given  by  the  product  of  the 
angular  error  in  local  hour  angle  (ALHA),  the  earth’s  radius  and  the  cosine  of  the  latitude. 
Using  Eq.  (F6) 


R(cos  L)  ALHA  = 


^cos  LHA  |/ (x?  +  Af^  sin  LHA^ 


(FIO) 


2  f^  o?  cos  8 


